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ABSTRACT 

Although point caustics harbour a larger potential for measuring the brightness profile of stars 
during the course of a microlensing event than (line-shaped) fold caustics, the effect of lens 
binarity significantly limits the achievable accuracy. Therefore, corresponding close-impact 
events make a less favourable case for limb-darkening measurements than those events that 
involve fold-caustic passages, from which precision measurements can easily and routinely 
be obtained. Examples involving later Bulge giants indicate that a ~ 10 % misestimate on the 
limb-darkening coefficient can result with the assumption of a single-lens model that looks 
acceptable, unless the precision of the photometric measurements is pushed below the 1 %- 
level even for these favourable targets. In contrast, measurement uncertainties on the proper 
motion between lens and source are dominated by the assessment of the angular radius of 
the source star and remain practically unaffected by lens binarity. Rather than judging the 
goodness-of-fit by means of a x 2 test only, run tests provide useful additional information 
that can lead to the rejection of models and the detection of lens binarity in close-impact 
microlensing events. 
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1 INTRODUCTION 

In order to resolve the surface of the observed star dur- 
ing a microlensing event, the magnification pattern cre- 
ated by the lens needs to supply a large magnification 
gradient. Two such configurations meeting this requirement 
have been discussed extensively in the literature: a point- 
causti c at the angular position of a single point-li k e lens 
<Witt & Mad 1 1994 iNemiroff & Wickramasinghef 1 1994 IGouTcI 
1994 iBoedanoy & Cherepashchukl Il995l 119961: IWitJ Il995t 
Gould &Welchl Il996t iGaudi & Gould Il999t iHevrovskv et all 
200G; Hevrovskv 2003) an d a line-shaped fold caustic pro duced 
by a b inary lens Ischn ridex&Weijjl^^ 

I l987t lGaudi& Gould Il999t iRhie & Bennettl 1 1 999t iDorninik 
2004a bid). It has been pointed out bv lGaudi & Gould <199Sl) that 
fold-caustic events are more common and their observation is eas- 
ier to plan, whereas close-impact events where the source transits a 
point caustic can provide more information. However, I will argue 
that this apparent gain of information can usually not be realized 
due to potential lens binarity. In contrast to fold caustics which form 
a generically stable singularity, point caustics are not stable and do 
not exist in reality. Instead, there is always a small diamond-shaped 
caustic containing four cusps. 

In this paper, the influence of lens binarity on the measure- 
ment of stellar limb-darkening coefficients and proper motion is in- 
vestigated and the arising limitations of the power of close-impact 
events where the source passes over a single closed caustic are dis- 



cussed. Sect. |2| discusses the basics of close-impact microlensing 
events with the effect of source size, the potential of measuring 
stellar proper motion and limb darkening, and the effect of lens bi- 
narity. Sect.[3]shows the influence of lens binarity on the extraction 
of information from such events. First, the effect of lens binarity on 
the light curves is demonstrated by means of two illustrative exam- 
ples involving K and M Bulge giants. Subsequently, a simulation of 
data corresponding to these configurations is used to investigate the 
potential misestimates of parameters if lens binarity is neglected. 
Sect.|4]presents the final conclusions and a summary of the results. 



2 CLOSE-IMPACT MICROLENSING EVENTS 
2.1 Size of source star 

As pointed out by IPaczvnskil ll986l) . a point-like source star at a 
distance Ds from the observer exhibits a magnification due to the 
gravitational field of a lens star with mass M at Dl by a factor 

a(u)= "!±jL , a) 

u yfv? + 4 

where source and lens are separated by the angle u 6e and 



4GM D s - D h 
denotes the angular Einstein radius. 
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The proper motion p of the source relative to the lens con- 
stitutes a microlensing event with the time-scale £e = #e/p, f° r 
which the lens-source separation becomes 



impact paramete r falls below the angul ar source radius, i.e. 
un < q+ iBosdanov & Cherepashchuld Il99d iGould & Welcf 



(3) 



u(t) = V«§ + b(*)] 2 , 

where p(t) = (t — to) /t^, so that no is the closest approach which 
occurs at time to- 

It was discussed bv Iwitt & Mad ll994l) as well as by 
iNemiroff & WickramasingheT^994l) that the finite extent of the 
source star could cause observable dev iations from t he magnifi- 
cation factor A(u) as given by Eq. fll lGould 1 1994) has argued 
that for a source star with radius R t — p* 6*e and radial brightness 
profile I(p) = I £(p), where p is the fractional radius and / is the 
average brightness, the finite-source magnification can be approxi- 
mated as 



A e *\u, Pt ;0 = A(u)B(u/p*;0- 



with 



2tt pi 

z _ 

y/p 2 + 2pz cos ip + . 



(4) 



-.p dp dip. (5) 



In particular, B(z) ~ 1 for z 3> 1, so that for large angular sepa- 
rations (compared to p*), the finite source effects becomes negligi- 
ble, whereas B(z) ~ /3z for z< 1 implies strong effects for small 
angular separations. Therefore, microlensing events with small im- 
pact parameters uo <C p are the most likely to show prominent 
effects of finite source size. Stellar spectra can be used to derive 
radius R* and distance Ds of the source star, yielding its angular 
radius 9* = R^/Ds- On the other hand, microlensing observations 
yield p* and therefore the time-scale t* = p* t%, during which the 
source moves by its own angular radius on the sky. Therefore, the 
observation of finite source effects in microlensing events provides 
a measurement of the proper motion between lens and source as 

R* 



^ U Ds p* t E 



2.2 Limb darkening 



(6) 



As already indicated, the stellar surface is not uniformly bright, 
but a characteristic variation with the distance from the center is 
observed, commonly known as limb darkening, which depends on 
wavelength and therefore on the filter used for the observati ons. A 
widely used model is the linear limb-darkening law lMilnell92ll) 



£(/>) = i + r (§vT 



-p 2 - 1 



(7) 



with ^ r ^ 1, which is linear in cosi? = y/l — p 2 , 
where i? is the emergent angle. If the point-source magnifica- 
tion A(u) shows a strong variation over the face of the source 
star, dense and precise microlensing observations provide an op- 
portu nity for a m easurem ent of t he limb-darkening coefficient 
teogdanov & Cherepashchukll996tlGould & Welchll99d) . 

The strongest magnification gradients occur in the vicin- 
ity of caustics. While a single lens creates a point caustic at 
its angular position, binary lenses create finite caustic curves 
which contain cusps. As the angular separation between the 
binary lens objects tends to zero, its diamond-shaped caustic 
with four cusps degenerates into the point caustic of a single 
lens. There are two different main scenarios with the potential 
of providing limb-darkening measurements: passages of the 
source star over the point caustic of a single lens, so that the 



1996} iGaudi & Gould Il999t IHevrovskv et alJ l200rA IHevrovskv 
2003), and passages of the source star ov er a (line-s haped) 
fold c austic created by a binary l ens {S chneider & WeiB 
[ 19871 ISchneider & Wagoned Il987t IGaudi & Gould Il999t 
iRhie & Bennett! Il999t iDominild l2004jd) . In addition, the 
source might pass di rectly over a cusp, a s for the event MA- 
CHO 1997-BLG-28 iAlbrow et all [T999b). for which the first 
limb-darkening measureme nt by m i crolens ing has been obtained. 
As anticipated by IGaudi & Gould ^999), the vast majority of 
other limb-darkening measurem ents so far has arisen from fold- 
caustic passages <Afonso et alj|200c| : lAlbrow et all l2000l 1200 it 
iFields etalj|2003l) . whereas two measurements from sing le-lens 
event s have been reported so far lAlcock et all 1 19971 lYoo et alJ 
120041) . The remaining limb-darkening measure ment by microlens- 
ing, on the solar-like star MOA 2002-BLG-33 lAbe et alJ <2003l) 
constituted a very special case, where the source simultaneously 
enclosed several cusps over the course of its passage. 



2.3 Lens binarity 

Since the majority of stars resides in some form of binary or mul- 
tiple systems (e.g. Abt 1983), it seems at first sight a bit surprising 
that more than 85 % of the observed microlensing events appear 
to be consistent with the assumption of a single point-like lens. An 
important clue to this puzzle is that the separation of binaries cov- 
ers a broad range of 6-7 orders of magnitude roughly from contact 
to typical stellar distances. Lens binarity however usually does not 
provide strong deviations to the light curve if the angular separation 
is much smaller or much la rger than the angular Einstein radius 
Be iMao & Paczvnsk3ll99ll) . Therefore, many binary lenses sim- 
ply escape our att ention by failing to provide an observable signal 
(Di Stefano 2000). While several binar y lens systems provide weak 
distor tions, such as MACHO LM C-1 iDominik & HirshfelJl994l 
1 1 9961) or MACHO 1999-BLG-47 jAlbrow et alj2002t) . a character- 
istic signature is provide d if the source passe s over a fold or even 
a cusp caustic, OGLE-7 1 Udals ki et alJI 19941) being the first such 
observed event. 

Thus, although a lens star is never a completely isolated ob- 
ject, it can be approximated as such in many cases. However, one 
needs to keep in mind that its caustic is a small diamond with 
four cusps rather than a single point. While for a single lens, the 
computation of the light curve is a rather straightforward process, 
where a semi-analytical expression b y means of e lliptical integral 
exists for uniformly bright sources I Witt & Mao 1994), it becomes 
quite demanding if a binary lens is considered due to the fact that 
no closed expression exists for the point-source magnification of 
a binar y lens, but in general a fifth-ord er polynomial needs to be 
solved iWitt & Madll995tlAsad3l2002h . For this paper, the algo- 
rithm of Dominik 1 1998) has been used which is based on the con- 
tour plot technique by Schramm & Kavser 1 1987) and the applica- 
tion of Green's theorem. 

Compared to a single lens, a binary lens involves three ad- 
ditional parameters, which can be chosen as the mass ratio q, the 
angle a, and the separation parameter d. Here, q = M2 jM\ is the 
ratio between the masses of secondary and the primary component 
which are separated on the sky by the angle d#E, while a is mea- 
sured from the vector pointing from the secondary to the primary 
towards the source trajectory which therefore reads 
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Table 1. Model parameters and indicative corresponding physical proper- 
ties for the two discussed event configurations. 





configuration I 


configuration II 


*E [d] 


55 


35 


*0 [d] 








itrj 


0.015 


0.05 


P* 


0.05 


0.075 


a 


o 




o 




9 


1 


1 


-^base 


13.6 


12.3 


9 








Ti 


0.5 


0.5 


D S [kpc] 


~ 8.5 


~ 8.5 


D L [kpc] 


~ 6.5 


~ 6.5 


M [M ] 
v [kms -1 ] 


~ 0.35 


~ 0.7 


~ 55 


~ 140 


u [uas d ] 


~ 330 


~ 460 


6» E [uas] 


~ 6 


~ 15 


r E [AU] 


~ 2.0 


~ 2.9 


^ [«©] 


~ 500 


~ 800 


R* [%] 


~ 25 


~ 60 


spectral type 


K5m 


M4III 


My 


-0.2 


+0.2 


V- I 


+2.0 


+3.4 


A r 


1.25 


0.85 


rn - M 


14.65 


14.65 



«(*) 



U0 



— sin a 

COSQ 



+ P(t) 



cos a 
sin a 



(8) 



3 BINARITY VS LIMB DARKENING 

3.1 Illustrative binary-lens configurations 

In order to investigate the effect of lens binarity on the measurement 
of stellar limb darkening and proper motion, two different event 
configurations have been chosen as illustrative examples, where 
different angular separations between the lens objects have been 
considered for both of them. For the model referred to as 'config- 
uration F, let us consider a K5 Bulge giant at Ds ~ 8.5 kpc, for 
which M v = -0.2 and V - I = +2.0. One therefore obtains 
lo = 12.35, so that with an assumed extinction Ai ~ 1.25, the ob- 
served source magnitude becomes / ~ 13.6. Let us further assume 
a binary lens of total mass M ~ 0.35 Mq at Dl ~ 6.5 kpc, so 
that the angular Einstein radius becomes 8e ~ 320 pas, the Ein- 
stein radius becomes te = Dl #e ~ 2.0 AU and its projection to 
the source distance becomes r' E = (Ds/Dl) te ~ 500 Rq. With 
Ri, ~ 25 Rq for a K5 giant, a source size parameter p* = 0.05 is 
therefore adopted. In order to make an optimal case for observing, 
a rather low proper motion fj, ~ 6 km s _1 , corresponding to a lens 
velocity of v ~ 55 kms^ 1 relative to the source has been chosen, 
which yields an event time-scale <e = 55 d. The impact parameter 
is chosen as uo = 0.015, corresponding to a peak magnification of 
a point source of Aq ~ 70. 

Configuration II has been chosen to be similar to the param- 
eters of the observed event MACHO 1995-BLG-30 lAlcock et all 
Il997t) . for which the source is an even larger star, namely an M4 
giant. According to the obtained model parameters, let us adopt 
uo = 0.05, p* = 0.075, and the event time-scale <e = 35 d. For 
Ds ~ 8.5 kpc and Dl ~ 6.5 kpc, the appropriate stellar radius 



of R* ~ 60 Rq is obtained for M ~ 0.7 Mq, so that r E ~ 
800 Rq. These choices yield r E ~ 2.9 AU and #e ~ 460 pas, 
so that the proper motion becomes p ~ 15 kms -1 and the rel- 
ative lens velocity is v ~ 140 kms" 1 . With My = +0.2 and 
V — I = +3.4 for an M4 giant, an extinction of Ai — 0.85 yields 
the baseline magnitude /base = 12.3. 

Binary lenses are likely to cause an asymmetry to the light 
curve, which however can be arbitrarily small and can even vanish 
for some configurations. In order to study the maximal impact of 
binarity on the measurement of limb darkening, configurations have 
been chosen that preserve the symmetry. Therefore, let us assume 
that both lens objects have the same mass and consider a source 
trajectory parallel to the line connecting their angular positions. 

In general, the light received from the source is blended with 
additional light from other unresolved sources (that are not affected 
by microlensing) or from the lens star, which is quantified by the 
blend ratio g — Fb/Fs, where Ds denotes the source flux and Db 
denotes the background (blend) flux. Since the choice of equal lens 
masses implies that the lens objects are M dwarfs of mass M/2 ~ 
0.18 M & or M/2 ~ 0.35 M Q , their contribution to the total light 
falls well below the systematic error bars even at the observed I- 
baseline. Therefore, blending is neglected with the choice g — 0. 
Finally, a limb-darkening coefficient 1+ = 0.5 has been adopted 
for both configurations. 

The model parameters and indicative corresponding physical 
properties of matching lens and source stars are summarized in Ta- 
ble □ Light curves that correspond to either of the adopted config- 
urations are shown in Fig. Q for a binary-lens separation parame- 
ter d — 0.2 along with light curves that correspond to single-lens 
models with otherwise identical parameters, while Fig.[5]shows the 
difference in magnitude between these lightcurves. Lens binarity 
decreases the magnification both around the peak and in the wing 
region, while an increase in magnification results in regions just be- 
fore the lens center is hit by the leading limb or just after it is hit 
by the trailing limb. For configuration I, these regions of increased 
magnification stretch over a much larger part of the light curve than 
the corresponding regions for configuration II. 

If we assume the lens center-of-mass as the origin of coordi- 
nates and the angular Einstein radius #e as the unitlength, and if 
we choose the separation vector of the binary lens objects along 
the a>axis, the four cusps of its central caustic are found at 



= (sj ,0) , 
s T = (As, st) , 



s^ = (s+ 0), 
Si = (As, -s 2 ) 



(9) 



For a point source, these four cusps degenerate into the single point 
s = 0, placing the caustic entry and exit at 



Po 



(*) = Va 



(10) 



The coordinates for which the leading or trailing limb of the 



source hits a cusp then follow from the condition [u(t) — s] 2 
with u(t) given by Eq. JSJ, as 



pt(t) 
Pt(t) 



pf(t) 



s 1 cos a ± \/ p\ — (uq + s 1 sin a) 



s~t cos a ± \J pi — (wo + s\ sin a 
S2 sin a + As cos a ± 



uq — S2 cos a + As sin a) 2 . 
S2 sin a + As cos a ± 
i-y/pf — (uq + S2 cos ol + As sin a) 2 . 



(11) 



For equal masses (q — 1), 



si and As = 0, while 
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t-f [d] t-t [d] 

Figure 1. Light curves corresponding to the binary-lens configuration I (left) and configuration II (right) with parameters as listed in TablefTland a separation 
parameter d = 0.2 (solid line) along with those for a single lens model with otherwise identical parameters (dashed line). In the lower panels, which show the 
full wing of the light curves, binary-lens and single-lens models are hardly distinguishable. 



for a point source, the four cusps degenerate into the point so = 0, 
placing the caustic entry and exit at 

Po(*) = Vp*- u o- (12) 

The passage of the leading stellar limb over the cusps is il- 
lustrated in Fig. [3] which shows snapshots of the source and the 
caustic for d = 0.2 and both configurations at times tZ-, tZ, tZ, or 
tZ , which correspond to the source positions p~ , p~, p7 , or p7 , 
respectively, along with the light curve for —4 d ^ t ^ 0. Dashed 
lines refer to tZ , corresponding to pZ , which is the beginning of 
the caustic passage for a single point lens. 

With the source size parameter being larger for configura- 
tion II, a smaller fraction of the source is subtended by the caustic 
for the same choice of the lens separation parameter d. A larger 
impact parameter no places the effect of binarity more towards the 
limb of the source. Despite a larger source size parameter p*, the 
combination of a smaller impact parameter uo and a smaller event 
time-scale t-E leads to a shorter caustic passage time, which, for a 
single lens, is tg — tZ ~ 3.9 d for configuration II compared to 
tg — tZ^ 5.2 d for configuration I. 

On the inclusion of a cusp, the leading limb may enter the 



caustic, exit the caustic, or just touch it there so that the tangents 
of the circle and the merging fold lines match. For d — 0.2 and 
configuration II, the caustic is entered for p~_ and p~ , whereas it is 
exited for p~ and pZ, . In contrast, the caustic is exited for pZ for 
configuration I and the same lens separation. A sign change in the 
curvature of the light curve similar to that for a single lens near t$ is 
seen for the binary-lens light curves near tjr, whereas all other cusp 
intersections are not that easily identified by observable features in 
the light curve. 

3.2 Simulated sampled light curves 

In order to see what kind of information can be extracted from sam- 
pled light curves, simulated data sets have been created that corre- 
spond to the two chosen binary-lens configurations. The event sam- 
pling parameters for this simulation have been chosen in analogy to 
an earlier investigation of fold-caustic passages iDomi nik 2004c). 
For data in the range t m i n ^ t ^ tmax, the sampling is character- 
ized by the sampling interval At, the fluctuation /a* of the time at 
which the measurement is taken, and the relative sampling phase 
shift /phase, where values of /At = 1/6 and /phase = 0.2 have 
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Figure 2. Difference between binary-lens light curves with the angular separation parameter d ■ 
their single-lens counterparts. 
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Figure 3. Passages of the leading limb of the source over the four cusps for configuration I (left) or configuration II (right) and a binary-lens separation 
of d = 0.2. The upper panels show the caustic (bold line) and the lens center (filled circle) with snapshots of the source as the trailing limb hits a cusp, 
as well as it touches the lens center (dashed line). The dotted line indicates the track of the source, while diagonal crosses mark the position when a cusp 
is hit, whereas the passage of the lens center is indicated by an open circle. For configuration I, the leading limb touches the caustic for p~ = —0.0672, 



P T = 
tz = 
for p~ 



-0.0497, p~ 



-0.0352, or pZ 



-0.0282 (in this order), corresponding to times 



-3.696 d, t~ 



-2.733 d, t 



-1.935 d, or 



-1.551 d relative to to, while the lens center is hit for p„ 
= -0.0754, p~ = -0.0690, p~ = -0.0364, orp^ 



-0.0477, i.e. att n 



-2.623 d. For configuration II, the leading limb touches the cusps 



to times t 



-2.639 d, t 



-2.4163 d, tZ 



-1.274 d, t 



-0.0255 (in this order), and hits the lens center for p ( y = —0.0559, which corresponds 
= —0.894 d, and t^ = —1.957 d relative to to. In contrast to configuration I, the order 
between and pZ and is reversed. The lower panels show the lightcurve for the binary lens (solid line) together with that for a point lens with otherwise 
identical parameters (dashed line), where t~ , t~ , t~ , and t~ are indicated by thin solid lines, and tjT is indicated by a thin dashed line. 
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Table 2. Sampling intervals and number of data points for configuration I 



tmin [d] 



[d] 



At 



N 






3 


20min 


405 


3 
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30 min 
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20 


lh 
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20 


30 


2 h 


226 


30 


60 


6h 


227 


60 


100 


12 h 


153 


fOO 


150 


Id 


95 


150 


300 


2d 


112 



With to = 0, a constant sampling rate has been adopted on both sides of 
the peak for intervals — t max ^ t ^ —tmin an d t m - ln ^ t ^ t max . Due 
to the applied phase shift /phase ar >d the fluctuation of the time when the 
observations are taken, /At, the actual times of observation may fall outside 
the originally designated interval. N denotes the number of data points for 
each selected t m ln and t max - In total, 2155 data points have been created. 



Table 3. Sampling intervals and number of data points for configuration II 
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Same as Tablel2lfor configuration II, where a total of 1741 data points have 
been created. 



been adopted here. The sampling intervals At for different regions 
of the light curve have been chosen to roughly comply with the ob- 
serving strategy of t he PLANET microlensing follow-up campaign 
I Albr ow et all 199 8: Dominik et al. 2002) for the central region and 
with that of the OGLE-III survey lUdalski 2003) for the outer re- 
gions. These are listed in Tables[2|and[3] 

For a reference magnitude of m ro f = 16 in J-band, an un- 
certainty of ovcf = 0.015 in the observed magnitude has been 
assumed, corresponding to a ~ 1.5 % relative uncertainty in the 
measured flux. For other magnitudes, it is assumed that the pho- 
tometric measurement follows Poisson statistics, so that its uncer- 
tainty is proportional to the square-root of the observed flux. This 
uncertainty has been smeared with a relative standard deviation of 
fcr = 0.125. The resulting photometric error bar <7 p hot y et does not 
represent the total measurement uncertainty. Instead, a systematic 
error errj is added in quadrature, which becomes dominant as the 
target gets bright. For the systematic error, three different values 
(To = 0.3 %, Co = 0.5 %, or (To = 1 % have been used. 

All the above choices mean that the creation of synthetic data 
sets(t (i) ,m (i) , 



a^) follows the rules 
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Figure 4. Binary-lens light curves (solid line) and simulated data for (To = 
1 % for both configurations and a lens separation parameter d = 0.2. Also 
shown are light curves for corresponding single-lens models with otherwise 
identical parameters (dashed line). The time tg at which the leading limb 
hits the single-lens point caustic is indicated by a vertical thin dashed line. 



N (m(t (i) ),cr (i) ) 



(17) 



where Af((i, cr) denotes a value drawn randomly from a normal 
distribution with mean fj, and standard deviation cr, m(v- 1 ') denotes 
the magnitude at time v-*' for the adopted model parameters, and 
i € [0, (t max train) /(At)]. 

In order to simulate possible losses in the data acquisition due 
to observing conditions or technical problems, data points have 
been removed from the dataset at random with a probability of 

Ploss = 5 %. 

In fact, a round-the-clock coverage (without daylight gaps) re- 
quires more than one telescope. Nevertheless, the effects of addi- 
tional free parameters (mainly baselines magnitudes) and different 
telescope characteristics for multi-site observations are neglected in 
favour of simplicity, and the data is treated such as resulting from a 
single telescope. 

The quality of coverage of the data set is demonstrated in 
Fig. |4] which shows the simulated data points and the theoretical 
light curves for d = 0.2 and ao — 1 % for both binary-lens config- 
urations around the contact of the leading limb with the lens center. 
To indicate the influence of lens binarity, the corresponding single- 
lens light curves are also shown. 

3.3 Best-matching single-lens models 

Let us now assume a single lens and determine the corresponding 
model parameters by means of fits to the simulated data sets that 
correspond to the two binary-lens configurations. We will then see 
whether statistical tests suggest to accept the single-lens model and 
how well the true limb-darkening coefficient V and the time-scale 
t*, which yields the proper motion as /j, — 6 t /t*, are reproduced. 
Any significant offsets for acceptable single-lens models will limit 
the accuracy to which these paramters are determined on the as- 
sumption of such a model. 

Tables|4]and|5]list the obtained model parameters and the re- 
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Table 4. Best-matching single-lens models to configuration I binary-lens models with different angular lens separations for different systematic errors. 







<t = 0.3 % 






<T = 0.5 % 






<T = 1 % 






d = 0.2 


d = 0.15 


d = 0.1 


d = 0.2 


d = 0.15 


d = 0.1 


d = 0.2 


d = 0.15 


d = 0.1 


X 2 ■ 

A min 


18054.9 


4337.6 


2235.6 


8873.9 


2969.9 


2097.89 


3895.1 


2267.8 


2037.6 


d.o.f. 


2149 


2149 


2149 


2149 


2149 


2149 


2149 


2149 


2149 




8.40 


2.02 


1.04 


4.13 


1.38 


0.976 


1.81 


1.06 


0.948 


\/ 2 xLn " V2n - 1 


124.5 


27.6 


1.32 


67.7 


11.5 


-0.78 


22.7 


1.80 


-1.72 


P(x 2 > xlJ 


< 10-3366 


< io- 167 


0.094 


< io-" 6 


< IO" 3 " 


0.781 


<io- 113 


0.037 


0.957 




55.44 


54.99 


54.93 


55.28 


54.94 


54.91 


55.18 


54.91 


54.91 


liO 


0.0157 


0.0149 


0.0148 


0.0155 


0.0148 


0.0147 


0.0153 


0.0148 


0.0147 


P* 


0.0541 


0.0518 


0.0503 


0.0542 


0.0519 


0.0503 


0.0543 


0.0520 


0.0503 


-^base 


13.609 


13.604 


13.601 


13.608 


13.604 


13.602 


13.607 


13.603 


13.602 


Yi 


1.000 


0.746 


0.538 


1.000 


0.751 


0.538 


1.000 


0.751 


0.537 



The complete data set has been included. The 'true' binary-lens parameters are shown in TablelTI The blending parameter has been fixed to g = 0, whereas 
to has been allowed to vary, yielding |<o| < 90 s in all cases. The quantity yj 2x 2 nin — \J2n — 1 is a measure of the goodness-of-fit as a characteristic of 
a x 2 -test yielding the equivalent deviation from the mean in units of the standard deviation for a Gaussian distribution, where n is the number of degrees 
of freedom (d.o.f.), which is the number of data points reduced by the number of free model parameters, while P(x 2 5 s Xmin) denotes the associated 
probability. 



Table 5. Best-matching single-lens models to configuration II binary-lens models with different angular lens separations for different systematic errors. 







(to = 0.3% 






(to = 0.5% 






(TO = 1 % 






d = 0.2 


d = 0.15 


d = 0.1 


d = 0.2 


d = 0.15 


d = 0.1 


d = 0.2 


d = 0.15 


d = 0.1 


x 2 ■ 


5670.3 


2371.2 


1791.5 


3161.4 


1895.0 


1666.8 


2001.6 


1668.7 


1605.0 


d.o.f. 


1735 


1735 


1735 


1735 


1735 


1735 


1735 


1735 


1735 


xlj<*.o.f. 


3.27 


1.37 


1.03 


1.82 


1.09 


0.961 


1.15 


0.962 


0.925 


V 2 xLn -V2n-1 


47.6 


10.0 


0.96 


20.6 


2.66 


-1.16 


4.4 


-1.13 


-2.24 


P(x 2 > xLJ 


< 10-494 


< IO" 23 


0.168 


< IO" 94 


0.004 


0.877 


< io- 5 


0.871 


0.988 




34.62 


34.76 


34.91 


34.62 


34.77 


34.91 


34.64 


34.78 


34.93 


"0 


0.0489 


0.0491 


0.0497 


0.0489 


0.0491 


0.0497 


0, 0489 


0, 0491 


0, 0497 


p* 


0.0817 


0.0795 


0.0758 


0.0817 


0.0796 


0.0758 


0.0816 


0.0796 


0.0757 


-^base 


12.302 


12.301 


12.301 


12.302 


12.301 


12.301 


12.302 


12.302 


12.301 




1.000 


0.984 


0.607 


1.000 


0.997 


0.605 


1.000 


1.000 


0.599 



The same quantities as in Tablel4lare displayed. Again, all data points have been included in the fits and the 'true' binary-lens parameters are listed in TablelTI 
The blending parameter has been fixed to g = 0, while to was allowed to vary, where | i o I i; 90 s resulted in all cases. 



suits of x tests for both binary-lens configurations and the differ- 
ent lens separations d and systematic errors uq. For configuration 
I, x tests recommend to accept the single-lens model for d = 0.1 
and all applied systematic errors, and marginally for d = 0.15 and 
a — 1 %. The limb-darkening coefficient TV turns out to be ~ 8 % 
larger for d = 0.1, ~ 50 % larger for d = 0.15, and tends to its 
maximal value Tr = 1 for d — 0.2. Compared to Yi, the rela- 
tive offsets on t+ are much smaller, namely ~ 0.5 % for d = 0.1, 
~ 4 % for d = 0.15, and ~ 10 % for d = 0.2, usually below the 
uncertainty in the measurement of . 

Similar to configuration I, the acceptance of the single-lens 
model for d = 0.1 and all applied systematic errors as well as for 
d — 0.15 and a — 1 % is recommended by \ 2 tests for config- 
uration II. A larger goodness-of-fit than for configuration I comes 
along with a larger offset in Yi, namely ~ 20 % for d = 0.1, while 
Yi ~ 1 is already reached for d = 0.15. The time-scale t* is found 
to be shifted by ~ 0.8 % for d = 0.1, ~ 5 % for d = 0.15, and 
~ 8 % for d = 0.2, 

The difference between the light curves for best-matching 
single-lens models and the 'true' binary-lens configuration for the 
three lens separations d — 0.2, d = 0.15, or d = 0.1 is shown in 



Fig.|U As Tables [4] and [5] show, the model parameters differ only 
very slightly for the different choices of the systematic error ao, 
so that the corresponding light curves are practically identical. For 
Fig-IU ao = 0.3 % has been chosen. For configuration I, some of 
the discrepancy is moved into the wing region of the light curve, 
where differences of up to 2 % (for d = 0.2), 0.7 % (for d = 0.15) 
or 0.2% (for d — 0.1) occur. In contrast, for configuration II, 
the maximal differences in the wing region are only 0.25% (for 
d = 0.2), 0.13% (for d = 0.15) or 0.1% (for d = 0.1). This 
reflects the fact that the region for which the binary-lens magnifi- 
cation exceeds that of a single lens extends significantly outside the 
region where the stellar limb touches the caustic for configuration 
I, while this is not the case for configuration II. The most signif- 
icant deviations occur around times when the stellar limb touches 
the caustic, which are however below 3 % (for d — 0.2), 2 % (for 
d — 0.15) or 1.2 % (for d = 0.1) for configuration I, and below 
2.5 % (for d = 0.2), 1 % (for d = 0.15) or 0.8 % (for d = 0.1) for 
configuration II. 

Since it is quite challenging to achieve photometric uncertain- 
ties below the 1 %-level, one should consider the effect of lens bina- 
rity as a serious problem for measuring limb-darkening coefficients 
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Figure 5. Differences in the light curve between a binary-lens model and a best-matching single lens model. The binary-lens models shown correspond to 
configuration I (left) or configuration II (right) and separations of d = 0.2, d = 0.15, and d = 0.1, while the single-lens model is the best fit for crn = 0.3 % 
as listed in TableRlor Tablel5l respectively. The upper panels show the full wing of the light curve, whereas the lower panels show only the peak region. 



from close-impact microlensing events and be quite careful with 
results obtained under the assumption of a single lens. 

An additional complication arises from the fact that photo- 
metric ^rrorbars reported by reduction algorithms such as DoPhot 
ISchechter et alll 9931 tend to underestimate the true photometric 
uncertainty for microlensing observations, which are plagued by 
varying observing conditions during the course of an event. Apart 
from the systematic error that is taken into account for the simu- 
lated data set, unde restimates of 10-20 % occur frequently (e.g. 
lUdalski et allll994l : lAlbrow et alj[l99ll iTsapras et alj|2003l) . so 
that the ratio x 2 /d.o.f. is not unlikely to exceed 1.4. This strongly 
limits the significance of the \ 2 test as a measure of the goodness- 
of-fit. 

If one accepts an underestimate by 20 %, the single-lens 
model for the configuration I binary-lens data with d — 0.15 and 
a = 0.5 % becomes acceptable, while for configuration II, single- 
lens models look acceptable for d = 0.15 and all applied system- 
atic errors as well as for d — 0.2 and a = 1 %. 

In addition to fits making use of the complete set of data 
points, best-matching single-lens models for the peak region of the 
lightcurve, defined as— 3 d ^ t ^ 3 d, have been obtained for 
(To = 0.3 %, for which the model parameters and the result of \ 2 
tests are displayed in Table|6| In contrast to the fits recognizing the 
full data set, the baseline /base has been fixed to its 'true' value, 
while the blending parameter g has been allowed to vary. As one 
could have expected from the larger discrepancies being attributed 
to the wing regions of the light curves for configuration I, the model 
parameters for the fits restricted to the peak region deviate more 
strongly from those obtained when the full data set is considered. In 
order to adjust to the optimal model, the blending parameter g has 
assumed a significant non-zero value. For all selected binary-lens 
separations d for configuration II and for d — 0.1 for configuration 
I, the goodness-of-fit resulting from the \ 2 test f° r the peak region 
model and data is worse than for the models and data for the full 
light curve, so that a rejection of all considered single-lens models 
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Figure 6. Differences in the light curve between a binary-lens model and 
a best-matching single lens model. Comparison between fits for the full 
data set (bold lines) and those where only the peak region (—3d ^ t ^ 
3 d) has been included (thin lines) with <jq = 0.3 % for both discussed 
configurations and a binary separation d = 0.1. 



is indicated. However, if one accounts for a possible 20 % increase 
in the size of the photometric errors, the models with d — 0.1 for 
both configurations still survive. In any case, the d = 0.1 models 
are not recommended for rejection for the larger systematic errors 
do = 0.5 % or (Jo = 1 % by means of a x 2 test over me peak 
region. 

Figure [6] shows a comparison of the differences between the 
binary-lens and the single-lens light curves of fits involving the full 
data set and those restricted to the peak region for a binary-lens 
separation parameter d = 0.1. While the differences between these 
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Table 6. Best-matching single-lens models to the peak region of binary-lens light curves for both discussed configurations and different angular lens separa- 
tions. 



<t = 0.3% 

configuration I configuration II 

d = 0.2 <2 = 0.15 d = 0.1 d = 0.2 d = 0.15 d = 0.1 



X, nin 1121.4 926.8 461.9 3468.21 948.98 572.34 

d.oT 398 398 398 398 398 398 

X^ in /d.o.f. 2.82 2.32 1.16 8.7 2.38 1.44 

\/ 2 Xmin ~ V 2n ~ 1 19 - 2 14 - 9 2 ' 20 55J 15 - 4 5 - 6 

PiX^X 1 ■ ) < 10- 81 < 10- 49 0.015 < 10~ 661 < 10- 53 < 10- 8 



*E 90.68 58.51 56.01 34.02 34.47 34.92 

u 0.0140 0.0181 0.0159 0.0478 0.0487 0.497 

p, 0.0332 0.0487 0.0495 0.0829 0.0802 0.0758 

g 0.571 0.031 0.011 0.0034 0.0014 0.0005 

Vj 1.000 0.663 0.531 1.000 1.000 0.605 



Results of fits of single-lens models to the peak region (— 3 d sC t sC 3 d) of the simulated binary-lens data sets for both configurations for a systematic error 
Co = 0.3 % and different angular lens separations. The parameters of the underlying binary-lens models can be found in Table HI In contrast to the fits that 
include the full data set, the baseline magnitude /base nas been fixed, whereas the blending parameter g has been allowed to vary. Again, the free parameter 
to fulfilled | to | < 90 s in all cases. 

has to be fulfilled. If one looks at the difference between single 
and binary-lens light curves, one sees that it becomes undetectable 
as soon as it is overshadowed by the statistical spread of the data 
which appears to be the case for d — 0.1 and ao = 1 %. In fact, 
the \ 2 test is n °t limited by this effect, but in principle allows to 
detect signals below the noise level by means of a sufficiently large 
number of independent measurements. 

The simulations show that single lens models for the discussed 
configurations involving a K or M Bulge giant look acceptable for 
d = 0.1, unless the photometric uncertainties are pushed signifi- 
cantly below the 1 %-level. The corresponding offset in the limb- 
darkening coefficient of ~ 10 % implies that the assumption of a 
single lens is incompatible with the desire of a precision measure- 
ment on r. In contrast, the effect on the time-scale t* is below 1 %, 
so that potential lens binarity does not have a significant effect on 
the measurement of the proper motion fi = 9* ft*, given that the 
uncertainty in is much larger. One might want to argue that the 
binary lens system is an unlikely configuration. However, d — 0.1 
means a projected separation a p = 0.2 AU for configuration I with 
two lens stars of mass M/2 ~ 0.18 Mq and a p = 0.3 AU for con- 
figuration II with two lens stars of mass M/2 ~ 0.35 Mq which 
should not a-priori be discarded. 



4 FINAL CONCLUSIONS AND SUMMARY 

Generally speaking, the potential binarity of the lens limits the ac- 
curacy of a limb-darkening measurement of the observed source 
star for close-impact microlensing events. The caustic near the cen- 
ter of the lens star is never exactly point-like, but always a small 
diamond with four cusps. If a point lens is assumed, the obtained 
limb-darkening coefficient is systematically offset, while the nature 
of the lens binarity may not be apparent. Although the inclusion of 
a binary lens in the modelling of the event will yield a proper limb- 
darkening measurement, the additional degrees of freedom still di- 
minish the achievable precision. Unfortunately, such a computation 
is extremely demanding. Simulations for typical event configura- 
tions involving K or M Bulge giants show that single-lens models 
for binary lens events that involve an offset of ~ 10 % on the limb- 
darkening coefficient T look acceptable both from x 2 and from run 



two types of fits are negligible for configuration II, the apparent dif- 
ferences for configuration I demonstrate the amount of information 
contained in the region before the first and after the last contact of 
the stellar limb with the caustic. For configuration I, the maximal 
deviation is reduced from 1.2 % to 0.8 % for the peak-only fits in 
the peak region at the cost of a deviation of up to 1.4% outside 
compared to 0.2 % when the full data set is considered. 

Given the problems in rejecting single-lens models by means 
of a x 2 test' it seems to be a good idea to look for characteristic de- 
viation patterns in order to decide on whether a model is in agree- 
ment with the observed data. An underestimate of error bars is not a 
problem for the assessment of a run test, which recognizes the sign 
of the residuals only, but not their size. In this sense, it is comple- 
mentary to the x 2 test, which is sensitive to the absolute value of 
the residuals, but blind to their signs. Although the obtained Xmin 
may look appropriate, a model is not acceptable if it fails a run test. 
Let a 'run' be defined as the longest contiguous sequence of resid- 
uals with the same sign, and let N denote the total number of data 
points, N+ the number of points with positive residuals, and Af- 
file number of points with negative residuals. For N > 10, the dis- 
tribution of the number of runs n r can be fairly approximated by a 
normal distribution with the expectation value 

£(n r )=l + ^ (18) 
and the standard deviation 



* ^ ^- if • (19) 

The goodness-of-fit can then be measured by the probability P T = 
P(n T ^ n° bs ). In addition to n T , other statistics such as the length 
of the longest run or the symmetry between N+ and jV_ may be 
checked. 

For the fits to the full data sets whose parameters are listed in 
Tables|4]and|5| Tables0and|8|show the result of the corresponding 
run tests over the full light curve and over the peak region only. Re- 
qiring an associated probabilitiy P T ^ 0.05 only lets the single-lens 
models for d = 0.1 and a — 1 % or a = 0.5 % for both configura- 
tions survive, while also the model for d = 0.15 and a = 1 % for 
configuation I nearly makes it to the acceptable region if P r ^ 0.01 
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Table 7. Results of run tests for the best-matching single-lens models for configuration I. 









a = 0.3% 






a = 0.5% 






(TO = 1 % 








d = 0.2 


d = 0.15 


d = 0.1 


d = 0.2 


d = 0.15 


d = 0.1 


d = 0.2 


d = 0.15 


d = 0.1 




N+ 


911 


1013 


1051 


935 


1034 


1063 


1016 


1053 


1071 




N- 


1244 


1142 


1104 


1220 


1121 


1092 


1139 


1102 


1084 


N = 2155 
(all data) 


S{fit) 
<r(n r ) 


1052.8 
22.7 


1074.6 
23.1 


1077.8 
23.2 


1059.7 
22.8 


1076.7 
23.2 


1078.3 
23.2 


1075.0 
23.1 


1077.9 
23.2 


1078.5 
23.2 


n° hs 


365 


774 


1064 


523 


928 


1086 


816 


1032 


1122 




S 


30.4 


13.0 


0.60 


23.5 


6.4 


-0.33 


11.2 


1.98 


-1.88 




1\ 


< IO" 202 


< IO" 38 


0.275 


< 1Q -122 


<io- 10 


0.630 


< IO" 28 


0.024 


0.970 




N+ 


187 


183 


197 


169 


179 


195 


164 


186 


199 




JV_ 


217 


221 


207 


235 


225 


209 


240 


218 


205 


N = 404 
(peak only) 


£{n T ) 
<r(n r ) 


201.9 
10.0 


201.2 
10.0 


202.9 
10.0 


197.6 
9.8 


200.4 
9.9 


202.8 
10.0 


195.9 
9.7 


201.7 
10.0 


203.0 
10.4 


n° bs 


111 


120 


190 


123 


152 


198 


152 


178 


204 




6 


9.1 


8.2 


1.28 


7.6 


4.9 


0.48 


4.5 


2.4 


-0.10 




Pt 


< io- 19 


< IO" 15 


0.100 


< io- 14 


< IO" 6 


0.317 


< 10~ 5 


0.009 


0.541 



For the fits of single-lens models to the full data set corresponding to configuration I and different values of the lens separation d and the systematic error <tq 
whose parameters are displayed in Table l4l run tests over the whole data set or the peak region (— 3 d $J t ^ 3 d) revealed JV+ positive and JV— negative 
residuals, where 'positive' means that the observed magnification exceeds the theoretical one. £ (n r ) denotes the expected number of runs and cr(n r ) denotes 
the corresponding standard deviation. For n° bs runs being found, the deviation in units of a(n r ) becomes 8 = [£(n r ) — ra° bs ]/cr(n r ), which corresponds 
to a probability P r = P(n r < n° bs ). 



Table 8. Results of run tests for the best-matching single-lens models for configuration II. 



CO = 0.3 % a = 0.5 % o-o = l % 







d = 0.2 


d = 0.15 


d = 0.1 


d = 0.2 


d = 0.15 


d = 0.1 


d = 0.2 


d = 0.15 


d = 0.1 




N+ 


888 


852 


846 


869 


850 


841 


856 


856 


839 




N— 


853 


889 


905 


872 


891 


900 


885 


885 


902 


N = 1741 
(all data) 


£(n v ) 
a{n r ) 


871.2 
20.9 


871.1 
20.9 


870.8 
20.8 


871.5 
20.9 


871.0 
20.9 


870.5 
20.8 


871.3 
20.9 


871.3 
20.9 


870.4 
20.8 


~ obs 
"r 


665 


776 


852 


729 


828 


884 


810 


862 


884 




S 


9.9 


4.6 


0.90 


6.8 


2.1 


-0.65 


2.9 


0.44 


-0.66 




Fr 


< 10~ 22 


<io- 5 


0.183 


< io- 11 


0.020 


0.742 


0.002 


0.329 


0.744 




N+ 


218 


189 


200 


203 


187 


200 


196 


198 


196 




N— 


186 


215 


204 


201 


217 


204 


208 


206 


208 


N = 404 
(peak only) 


£(n T ) 
<x(n r ) 


201.7 
10.0 


202.2 
10.0 


203.0 
10.0 


203.0 
10.0 


201.9 
10.0 


203.0 
10.0 


202.8 
10.0 


202.9 
10.0 


202.8 
10.0 


~obs 
"r 


51 


107 


166 


93 


147 


190 


141 


174 


194 




S 


15.1 


9.5 


3.7 


11.0 


5.5 


1.29 


6.2 


2.9 


0.88 




Fr 


< 10~ 51 


< IO" 21 


< 10" 4 


< io- 27 


< io- 7 


0.098 


< nr 9 


0.002 


0.189 



Results of run tests corresponding to fits of single-lens models to the full data set corresponding to configuration II and different values of the lens separation 
d and the systematic error <ro whose parameters are displayed in Tablel5l where either the whole data set or the peak region (— 3 d ^ t ^ 3 d) has been used. 
The quantities shown are the same as in TablelTl 



tests, unless photometry significantly below the 1 %-level is possi- 
ble. The measurement uncertainty of the proper motion /i = 8+/t+ 
however is dominated by the determination of the angular source 
radius involving stellar spectra, compared to which the accuracy 
limits on t* caused by lens binarity for close-impact microlensing 
events are negligible. 

In contrast, the measurement of limb darkening from fold- 
caustic passages does not suffer from any of the problems encoun- 
tered for close-impact events. A binary lens is assumed a-priori, 
but the measurement of limb darkening only depends on a smaller 
set of local properties rather than on the complete binary-lens pa- 
rameter space. The use of the local approximatio n of the light 
curve i n the vicinity of the fold-caustic passage (e.g. Albr ow et alJ 
Il999d : lDominikll 2004b ) makes the computation quite inexpensive 



and easy. Moreover, fold-caustic exits can be well predicted. Once 
a caustic entry has been observed, it is clear that a corresponding 
exit will occur, and the light curve on the rise to the caustic exit 
peak allows to predict the time of the caustic exit with sufficient 
precision usually more than a day in advance. A fair coverage of 
the caustic entry as well as the determination of the spectral type 
of the source star and an early measurement of the event time-scale 
tE even allows a rough guess on the passage duration and there- 
fore a proper a-priori assessment of the potential for measurements 
of the source brightness profile. While the characteristic proper- 
ties of a fold-caustic exit, which provides a full scan of the source 
from the leading to the trailing limb, can therefore be estimated 
in advance, corresponding predictions for close-impact events are 
only possible after the leading limb has alreaded passed the caus- 
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tic, leaving only the second half of the caustic passage involving the 
trailing limb. From fold-caustic passages, a linear limb-darkening 
coefficient can be obtained routinely and easily with a precision of 
less than 5 % (e.g. Dominik 2004c). Hence, events involving fold- 
caustic passages turn out to be the clear favourite over close-impact 
events for measuring limb-darkening coefficients apart from those 
where the source transits a cusp. The latter provide an opportunity 
to determine further limb-darkening coefficients beyond a linear 
law, which is quite challenging and usually impossible for fold- 
passage events lDominikl2004al) . 

However, events with fold-caustic passages make the worse 
case for proper motion measurements. Fits to the light curve in the 
vicinit y of the caustic passage only yield tjr = t*/ (si n <j>) (e.g. 
lAfonso et alJl998U Albrow et alll999c1 : lDominikl2004bl) . where <p 
is the caustic-crossing angle, which needs to be determined from 
a model involving the full set of binary lens parameters. While 
a precise measurement of t* is possible, the uncertainty in t* is 
severely limited by degener acies and appararent ambiguities for 
the caustic-crossin g angle <f> (Dominik 1999; Al brow et all 1999a: 
lAfonso et al1l200Cl) . which can have a comparable or even larger 
influence on the determination of the proper motion /i than uncer- 
tainties in the angular stellar radius 9*. 
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